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Abstract 

\ We compute the leading, chiral charge-changing relaxation term in the quantum transport 

' equations that govern electroweak baryogenesis using the closed time path formulation of non- 

CN ■ equilibrium quantum field theory. We show that the relaxation transport coefficients may be 

O . resonantly enhanced under appropriate conditions on electroweak model parameters and that such 

Q , enhancements can mitigate the impact of similar enhancements in the CP- violating source terms. 

\ We also develop a power counting in the time and energy scales entering electroweak baryogenesis 

' and include effects through second order in ratios e of the small and large scales. We illustrate 

the implications of the resonantly enhanced 0{e^) terms using the Minimal Supersymmetric Stan- 
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dard Model, focusing on the interplay between the requirements of baryogenesis and constraints 
Tij" , obtained from collider studies, precision electroweak data, and electric dipole moment searches. 
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I. INTRODUCTION 



The origin of the baryon asymmetry of the Universe (BAU) remains an important, un- 
solved problem for particle physics and cosmology. Assuming that the Universe was matter- 
antimatter symmetric at its birth, it is reasonable to suppose that interactions involving ele- 
mentary particles generated the BAU during subsequent cosmological evolution. As noted by 
Sakharov fi\ , obtaining a nonzero BAU requires both a departure from thermal equilibrium 
as well as the breakdown of various discrete symmetries: baryon number (B) conservation, 
charge conjugation (C) invariance, and invariance under the combined C and parity (P) 
transformations^. The Standard Model (SM) of strong and electroweak interactions satisfies 
these conditions and could, in principle, explain the observed size of the BAU: 



where is the baryon number density, s is the entropy density of the universe, and where 
the values shown correspond to 95% confidence level results obtained from Big Bang Nucle- 
osynthesis (BBN) and the cosmic microwave background (WMAP), respectively. In practice, 
however, neither the strength of the first-order electroweak phase transition in the SM nor 
the magnitude of SM CP-violating interactions are sufficient to prevent washout of any net 
baryon number created by P-violating electroweak sphaleron transitions during the phase 
transition. 

The search for physics beyond the SM is motivated, in part, by the desire to find new 
particles whose interactions could overcome the failure of the SM to explain the BAU. From 
a phenomenological standpoint, a particularly attractive possibility is that masses of such 
particles are not too different from weak scale and that their interactions both strengthen 
the first-order electroweak phase transition and provide the requisite level of CP- violation 
needed for the BAU. Precision electroweak measurements as well as direct searches for new 
particles at the Tevatron and Large Hadron Collider may test this possibility, and experiment 
already provides rather stringent constraints on some of the most widely considered exten- 
sions of the SM. In the Minimal Supersymmetric Standard Model (MSSM), for example, 
present lower bounds on the mass of the lightest Higgs boson leave open only a small window 
for a sufficiently strong first-order phase transition, althourfi this constraint may be relaxed 
by introducing new gauge degrees of freedom (see, e.g., y,|5|). Similarly, limits on the per- 
manent electric dipole moments (EDMs) of elementary particles and atoms imply that the 
CP- violating phases in the MSSM must be unnaturally small (~10~^). Whether such small 
phases (supersymmetric or otherwise) can provide for successful electroweak baryogenesis 
(EWB) has been an important consideration in past studies of this problem. 

In order to confront phenomenological constraints on the parameters of various elec- 
troweak models with the requirements of EWB, one must describe the microscopic dynamics 
of the electroweak phase transition in a realistic way. Theoretically, the basic mechanism 
driving baryogenesis during the phase transition is well-established. Weak sphaleron tran- 
sitions that conserve B — L but change B and L individually are unsuppressed in regions 
of spacetime where electroweak symmetry is unbroken, while they become exponentially 



Allowing for a breakdown of CPT invariance relaxes the requirement of departure from thermal equilib- 
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suppressed in regions of broken symmetry. Net baryon number is captured by expanding 
regions of broken symmetry ("bubbles"). Given sufficiently strong C and CP- violation as 
well as departure from thermal equilibrium, the non-zero B generated outside the bubble 
cannot be entirely washed out by elementary particle interactions that occur at the phase 
boundary. The baryon number density, Pb, is governed by a diffusion equation of the form: 

dtpB{x) - DV'^pb{x) = -T„sF„s{x)[nL{x) + Rpb{x)] , (2) 

where D is the diffusion coefficient for baryon number, r„s is the weak sphaleron transition 
rate, F^^six) is a sphaleron transition profile function that goes to zero inside the regions 
of broken electroweak symmetry and asymptotically to unity outside, i? is a relaxation 
coefficient for the decay of baryon number through weak sphaleron transitions, and nL{x) 
is the number density of left-handed doublet fields created by "fast" chirality changing 
processes (see, e.g., p). Thus, in order to obtain nonzero pb inside the bubble of broken 
electroweak symmetry, the left-handed density ni must be non- vanishing in the plasma at 
the phase boundary and possibly beyond into the region of unbroken symmetry. 

In effect, ni{x) acts as a seed for the S-changing weak sphaleron transitions, and its 
spacetime profile is determined by the CP- violating sources and the quantum transport of 
various charges in the non-equilibrium environment of the plasma. Typical treatments of 
these dynamics involve writing down a set of coupled quantum transport equations (QTEs) 
for the relevant charges, estimating (or parameterizing) the relevant transport coefficients, 
and solving the system of equations under the appropriate boundary conditions. 

Among the developments in the past decade or so which have made significant impacts 
on this program, we identify two that form the basis of our investigation in this work. 
First, the authors of Ref. [7| noted that diffusion of chiral charge ahead of the advancing 
phase transition boundary into the region of unbroken symmetry could enhance impact of 
baryon number-changing sphaleron processes, thereby leading to more effective EWB. The 
second, perhaps less widely-appreciated, development has been the observation by the author 
of Ref. ji| that the application of equilibrium quantum field theory (QFT) to transport 
properties in the plasma is not necessarily appropriate. In contrast to equilibrium quantum 
dynamics, the time evolution of quantum states during the phase transition is non-adiabatic. 
Consequently, scattering processes that drive quantum transport are no longer Markovian, 
but rather retain some memory of the system's quantum evolution. Using the closed time 
path (CTP) formulation of non-equilibrium QFT 9] to cornpute the CP-violating source 
terms in the plasma for the MSSM, the author of Ref. [8| found that these "memory 
effects" may lead to significant resonant enhancements (of order 10^) of the sources over 
their strength estimated in previous treatments (see, e.g., Ref. and references therein). 
The authors of Ref. ll|, subsequently found that performing an all-orders summation 



of scattering from Higgs backgrounds reduces the size of the CP- violating sources to some 
extent, but that the resonant enhancements nonetheless persist. Taken at face value, these 
enhancements would imply that successful EWB could occur with significantly smaller CP- 
violating phases than previously believed, thereby evading the present and prospective limits 
obtained from EDMs. 

To determine whether or not such conclusions are warranted, however, requires that one 
treat the other terms in the transport equations in the same manner as the CP- violating 
sources. Here, we attempt to do so, focusing on the terms that, in previous studies, have 
governed the relaxation of ni(x). In particular, chirality- changing Yukawa interactions 
with the Higgs fields and their spacetime varying vacuum expectation values (vevs) tend to 
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wash out excess nL{x). In earlier studies — including those in which non-equilibrium QFT 
has been applied to the CP- violating sources — these relaxation terms were estimated using 
conventional quantum transport theory 11, 3|- However, if the memory effects that 



enhance the CP- violating sources have a similar effect on these Yukawa terms, then the net 
effect on pb may not be as substantial as sugg ested in Refs. BlUQ- 

The goal of the present study is to address this question by developing a more compre- 
hensive treatment of EWB using the CTP formulation of non-equilibrium QFT. In doing 
so, we follow the direction suggested in Ref. Q and compute the transport coefficients of 
the chiral charges using the CTP formalism. To make the calculation more systematic, we 
identify the relevant energy and time scales that govern finite temperature, non-equilibrium 
dynamics and develop a power counting in the ratios of small to large scales (generically 
denoted here as e). As we show below, both the CP- violating sources and the driving relax- 
ation terms first arise at O(e^), and we truncate our analysis at this order. In contrast to the 
computation of the CP- violating sources, the derivation of the relaxation terms requires the 
use of finite density Green's functions. Given the resulting complexity, we consider here only 
the terms in the transport equations that previous authors have considered the dominant 
ones, and use our analysis of these terms to illustrate a method for obtaining a more com- 
prehensive treatment of the QTEs. To make the phenomenological implications concrete, 
we focus on the MSSM, realizing, however, that one may need to include extensions of the 
MSSM in order to satisfy the requirements of a strong first-order phase transition. Finally, 
we also attempt to identify the different approximations that have entered previous treat- 
ments of EWB, such as the implicit truncation at a given order in e and outline additional 
calculations needed to obtain a comprehensive treatment. 

Based on our analysis, we find that under that same conditions that lead to resonant 
enhancements of the CP- violating sources, one also obtains a similar, resonant en- 

hancement of the driving chirality-changing transport coefficient, F. Since Yb ^ S'^/Vt, 
resonant relaxation counteracts the enhanced sources, though some overall enhancement of 
EWB still persists. Consequently, it will be important in future work to study the other 
transport coefficients whose impact has been considered sub-leading, since they may be 
enhanced under other conditions than for the leading terms. From the standpoint of phe- 
nomenology, we also illustrate how the implications of EDM searches for EWB depends in a 
detailed way on the electroweak model of interest as well as results from collider experiments 
and precision electroweak data. 

In presenting our study, we attempt to be somewhat pedagogical, since the methods are, 
perhaps, not generally familiar to either the practitioners of field theory or experimentalists. 
Most of the formal development appears in Sections ITTVllVl In Section IH] we review the 
CTP formalism and its application to the QTEs and discuss in detail the formulation of 
density-dependent Green's functions. In Section IIIII we compute the CP- violating source 
terms, providing a check of Ref. P|, as well as the transport coefficients of the chiral charge 
densities. Here, we also enumerate the approximations used to obtain a set of coupled, linear 
differential diffusion equations, discuss their limits of validity, and identify additional terms 
(usually assumed to be sub- leading) that we defer to a future study. In Section |TV| we solve 
these equations for the baryon density. A reader primarily interested in the phenomenological 
implications may want to turn directly to Section |3 which gives illustrative numerical 
studies using the parameters of the MSSM. A discussion of the implications for EDMs also 
appears here. Section IVll contains a summary and outlook, while several technical details 
appear in the Appendices. 
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II. NON-EQUILIBRIUM TRANSPORT: CTP FORMULATION 



In what follows, we treat all CP-violating and non-topological chirality-changing inter- 
actions perturbatively^. In contrast to zero-temperature, equilibrium perturbation theory, 
however, the perturbative expansion under non-equilibrium, T > conditions requires the 
use of a more general set of Green's functions that take into account the non-adiabatic evo- 
lution of states as well as the presence of degeneracies in the thermal bath. Specifically, the 
matrix element of any operator 0{x) in the interaction representation is given by: 

{n\SlT{0{x)S,^,}\n) , (3) 

where 

Sint = T exp I z / d'^x Cint I (4) 



for an interaction Lagrangian C-mt, T is the time-ordering operator, and \n) is an in-state. 
In ordinary, zero-temperature equilibrium field theory, the assumptions of adiabaticity and 



of non- degeneracy of the states \n) implies that the only impact of 5*;^^ is the introduction 
of an overall phase, allowing one to rewrite Q as: 

{n\T{0{x)S,,,t}\n) 



{n\Sint\n) 



(5) 



This simplification is no longer valid for non-equilibrium T > evolution, and one must take 
into account the action of Sj^^ appearing to the left of 0{x) in Q. Doing so is facilitated by 
giving every field in S'int and Sj^^ a and "— " subscript respectively. The matrix element 
in then becomes: 

{n\V exp i (^j d^x - j '^-^ | 1^) . (6) 

where the path ordering operator V indicates that all "+" fields appear to the right of all 
"— " fields, with the former being ordered according to the usual time-ordering prescription 
and the latter being anti-time-ordered [here, 0{x) has been taken to be a "+" field]. Note 
that the two integrals in the exponential in © can be written as a single integral along a 
closed time path running from — oo to +oo and then back to — oo. 

Perturbation theory now proceeds from the matrix element © along the same lines as 
in ordinary field theory via the application of Wick's theorem, but with the more general 
V operator replacing the T operator. As a result, one now has a set of four two-point 
functions, corresponding to the different combinations of "+" and "— " fields that arise from 
contractions. It is convenient to write them as a matrix G{x,y): 

G(xv)- ( ^'("^'^^ -G<{x,y)\ , . 

lG>(x,y) -G\x,y) J 



^ Sphaleron transitions, however, are manifestly non-perturbative, and we parameterize their effects in the 
standard way. 
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where 



G>(x,i/) = (0_(x)0Uy)) 

G\x, y) = {T{<P+{x)<P+{y)]) = O^x^ - yo)G>{x, y) + e{y^ - x^)G<{x, y) 
G\x,y) = {T{(l>_{xW_{y)]) = e{x^~y^)G<{x,y) + 9{y,-x,)G>{x,y) 



(8a) 
(8b) 
(8c) 
(8d) 



and where the ( ) denote ensemble averages, 



(0(x)) = iTr[pO(x)] . 



(9) 



p is the density matrix containing information about the state of the system. In thermal 
equilibrium p is time-independent and is given by p = e"^'^'^"'^''^') for a grand-canonical 
ensemble. Note that the matrix G{x,y) may be written more compactly as: 



The presence of the factor is a bookkeeping device to keep track of the relative minus sign 
between the £+ and £_ terms in Eq. (jH)). 

The path-ordered two-point functions satisfy the Schwinger-Dyson equations: 



where the "0" superscript indicates a non-interacting Green's function and where E(x, y) is 
the matrix of interacting self energies defined analogously to the G{x,y). An analogous set 
of expressions apply for fermion Green's functions, with an appropriate insertion of —1 to 
account for anticommutation relations. 

A. Quantum Transport Equations from CTP Formalism 

The Schwinger-Dyson Eqs. (jlip are the starting point for obtaining the transport equa- 
tions governing ^^(x). To do so, we follow Ref. and apply the Klein-Gordon operator to 
G{x,y). Using 




(10) 




(lib) 



(11a) 



(a + m') G'{x, y) = [Uy + m") G\x, y) = {x - y) 



(12) 



gives 




(13b) 



(13a) 
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It is useful now to consider the (a, 6) = (1, 2) components of these equations: 



(a. + m2) G<{x,y) 



T.\x,z)G<{z,y)-T.<{x,z)G\z,y) 



{Uy + m^)G<{x,y) = -i J d^z ^G\x, z)i:<{z,y) - G<{x, z)J:\z,y) 

Subtracting Eq. ()14b|) from Eq. p4a|) and muhiplying through by i gives 
I - G<ix, y) = I hm 9j - d^) G<(x, y) , 

where X = {x + y) /2. However, 



\\m{cr- dl)G<{x,y) = -ij,{X) 



x->y 



(14a) 
(14b) 

(15) 

(16) 



where jfj.{x) = z(: 0^(a;)9^0(a;) :) = {n{x) , j (x)) , since the "+" and "— " labels simply indicate 
the order in which the fields (f)^{y) and (f){x) occur and may be dropped at this point. Finally, 
expressing G*'*(x, y) and S*'*(a:, ?/) in terms of ^-functions as in Eqs. ((HI), we obtain from 
Eq. (US): 



dn 



+ V-j(X) 



d^z 



dzo J:>{X,z)G<{z,X) -G>{X,z)i:<{z,X) 

+G<(x, z)j:>{z,x) -j:<{x,z)g>{z,x) 



(17) 



Following similar steps, but taking the sum rather than the difference of the components of 
the Schwinger-Dyson equations involving the S''{x,y) component on the LHS, one obtains 
the analogous continuity equation for Dirac fermions: 



dn 
dX'o 



V-j(X) = - dh 



dzn Tr 



S>(X, z)S<{z,X)- S>{X, ;z)S<(2,X) 



+S<{X,z)T>{z,X) - T.<{X, z)S>{z,X) 



(18) 



where 



Sai3{x^y) = {i^~c{x)i)+p{y)) (19a) 

Sap{x^y) = -{^-i3{y)^+a{x)) , (19b) 

displaying explicitly the spinor indices a, (3. Note that the overall sign of the RHS of Eqs. (jl7t 
ITH|) differs from that in Ref. ^ since the definition of our Green's functions G{x,y) and 
S{x,y) differ by an overall factor of —i. 

In many extensions of the SM, one encounters both chiral and Majorana fermions, which 
carry no conserved charge. It is useful, therefore, to derive the analogous continuity equation 
for the axial current jfj,5{x) = {ip{x)'jfj,j5ip{x))- Doing so involves multiplying the Schwinger- 
Dyson equations by 75, performing the trace, and taking the difference rather than the sum 
of the components involving S^{x,y) on the LHS. The result is: 

^ + V-j5(X) =2imP{X) (20) 



dX. 



+ / d^z 



Xo 



dzn Tr 



S>(X,2)5<(z,X) +5>(X, z)S<(2,X) 
- S< (X, z) S> (^, X) - S< (X, z) S> (^, X)] 75 } 
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where P{x) = {ip{x)'~f5ip{x)) and m is the fermion mass. In principle, one could evaluate 
P{x) using path-ordered perturbation theory as outlined above. 

B. Power Counting of Physical Scales 

Evaluating the various terms in Eqs. (fTTl ITHj) leads to a system of coupled quantum 
transport equations for the charges that ultimately determine ^^(x). On the LHS of these 
equations, it is conventional to parameterize j = — i5(Vn), in terms of the diffusion coefficient 
D (whose expressions we take from Ref. ^3). The RHS involves integrating the products of 
various Green's functions and self-energies over the entire history of the system. In practice, 
this integral depends on the various time and energy scales associated with non-equilibrium 
dynamics at finite temperature and density. Here, we observe that there exists a hierarchy 
among these scales that leads to a natural power counting in their ratios (generically denoted 
here as e) and that provides for a systematic expansion of the RHS of the transport equations 

(HIllliEnD. 

The changing geometry associated with the expanding region of broken symmetry and 
the spacetime variation of the Higgs vevs leads to a decoherence of states that have, initially, 
precise energy and momentum. The effect is analogous to the quantum mechanical evolution 
of a particle in a box of side L. If the value of L is changed to L+AL in some time interval At, 
a state that is initially a stationary state for the original box will become an admixture of the 
stationary states of new box. The shorter the interval At or the greater the wavenumber k 
of the initial state, the smaller the probability will be of finding the particle in the state with 
the same wavenumber in the new system. The time scale that characterizes this decoherence, 
Td, is naturally given by ~ 1/vk, where v = AL/At is the velocity of expansion of the 
box and k = p/h. In the present case, the relevant velocity is just v^, the expanding bubble 
wall velocity, the relevant wavenumber is = |k|. The smaller the velocity or the longer 
the wavelength, the more adiabatic the dynamics of the expanding bubble become and the 
longer the decoherence time. Equilibrium dynamics are approached in the adiabatic limit: 
Td ^ oo. The need to employ the OTP formalism follows from being in a situation with 
Vw > 0, or Td < oo. 

A second time scale that one encounters in quantum transport at the phase boundary 
arises from the presence of degeneracies among states in the thermal bath that vanish in 
the T — > limit. At finite T, for example, a single, on-shell fermion may be degenerate 
with another state involving an on-shell fermion-gluon pair — a situation that is forbidden at 
T = 0. Interactions of strength g that cause mixing between such degenerate states give rise 
to thermal — or plasma — widths Tp of order aT with a = g^/47r, and transitions between 
the degenerate states take place on a plasma time scale Tp of order ~ l/Ep. Again, the use 
of the OTP formalism is necessitated when Tp < oo or T > 0. 

A third time scale, which we denote Tint, is associated with the intrinsic frequency Uk of the 
quasiparticle states that characterize the plasma dynamics. This time scale is naturally given 
by Tint ~ l/ci^fc- In the present case, we note that although the decoherence and plasma times 
are finite, they are typically much smaller than Tint. For example, Tint/Td = v^k/uk < Vw/c. 
Numerical studies indicate that Vw/c <^ 1. Similarly, rint/rp = aT/uk- Since quasiparticle 
thermal masses are of order gT or larger, one also has that the latter ratio is smaller than 
unity. Thus, one is naturally led to expand the RHS of the transport equations in these 
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ratios: 



< Tint/Td < 1 (21a) 

< Ti^t/Tp < 1 . (21b) 

Finally, we observe that the generation of baryon number takes place in an environment 
of finite, but small particle number (or chiral charge) densities that are associated with 
chemical potentials /i,. For the temperatures and densities of interest here, one has |/Xi|/T <^ 
1, so that the latter ratio also provides for a natural expansion parameter. Denoting each 
of the ratios^ in Eq. ()2H) and /ij/T by e, we show below that both the CP-violating sources 
and the relaxation term first arise at (9(6^), and we truncate our analysis at this order. We 
note that doing so introduces some simplifications into the evaluation of the RHS of the 
transport equations. For example, both the self energies and the Green's functions G^, 
etc. depend on thermal distribution functions /(T, /ij) that differ, in general, from their 
equilibrium values, /o(T, /ij). The difference 6f = f(T,fii) — /o(T, /Xj) that characterizes 
the departure from equilibrium will be at least of 0{e), since it must vanish in the f u, — > 
limit. We find below that the effect of having 6f ^ contributes at higher order in e than 
we consider here, so that we may use the equilibrium distribution functions in the Green's 
functions and self-energies. 



C. Green's Functions at Nonzero Temperature and Density 

The computation of the various components of G{x, y) and S(x, y) appearing in Eqs. (fT7| 
ITHj) at nonzero temperature and density requires knowledge of (T, /ii)-dependent fermion and 
boson propagators. The T-dependence of propagators has been studied extensively (see, e.g. 
Ref. |1^J and references therein), while the /ij-dependence of fermion propagators has been 
studied in Refs. Here we summarize the features of (T, /ij)-dependent propagators which 
are important for our subsequent application of the real-time, CTP formalism of Sec. IHl and 
give some more technical details in Appendix ^ 

For pedagogical purposes, we provide here a brief derivation of the non-interacting fermion 
propagator but only give final results for the case of interacting fermions and bosons. To 
do so, we start from the mode expansions for the field operators appearing in the free Dirac 
Lagrangian, ipi^x) and %l){x): 

E [^«(k)«(k, a)e-^'- + 4(k)t;(k, a)^'-] (22a) 
E K(k)«(k, «)e^'^- + rf„(k)t;(k, a)e-^'-] , (22b) 



where = (cjk, k), uj^ = \/ |k|^ + m^, the mode operators satisfy: 



{6„(k),6T(k')} = {rf„(k),4(k')} = (27r)35(3)(k_k')2cuk5,^, (23) 



For our purposes, it is not necessary to distinguish a hierarchy among the different scale ratios. 
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and 



(24a) 
(24b) 



with f{uj,fii) being the non-equihbrium Fermi distribution function. For our purposes, 
the relative change 6 f {u , fii) / fo{u , fii) enters the transport equations multiplying exphcit 
factors of Tp and either or fi, so that in working to second order in e we may replace / 



by the equilibrium distributions /o(co',yUj 



+ 1] ^ Using the mode 



expansion it is straightforward to show that S^{x,y) = {ilj{x)tp{y)) and S^{x,y) 
— {ip{x)ijj{y)) can be expressed as (A denotes either ">" or "<"): 



S\x,y) 



d'^k 
(2^ 



-ik-{x—y) X 



gF{ko,IJ'i)p{ko,k) 



m] 



in terms of the free particle spectral density: 



p(^o,k) 



2uk 

and the functions: 



1 



ko — 'ji!k + i€ ko + Uk + it 



1 



ko — Uk — it fco + cu/; — it 



9F{ko,p-i 
gpiko^P-i 



1 - np^kQ - fii) 
-npiko - fii) . 



(25) 



(26) 



(27a) 
(27b) 



The propagators S'*'*(x, y) can now be constructed from the S''^(x, y) as in Eqs. (jH)). 

In the presence of interactions (characterized by a generic coupling g), the fermion prop- 
agator becomes considerably more complicated than given by Eq. ()25p. In particular, single 
fermion states can mix with other multiparticle states in the thermal bath, leading to the 
presence of additional poles (the "hole" modes) in the fermion propagator 0, 0|. The 
general structure of the fermion propagator arising from these effects has been studied ex- 
tensively at zero density In Appendix A we generalize to the case of non-zero For 
massless fermions, the resulting propagators are given by: 



S^{x,y;iJi 



d^k 
f27H 



-ik-{x—y) A 



9Fi.ko,fj) 



7o-7-k , . 7o + 7-k 
p+(fco, k, fii) H 



p_(A;o,k,/ii 



where k is the unit vector in the k direction, and 



p+(/co,k,yUi 



ko £p{k, fii 

Zhjk, —fii] 

ko + Sh{k,-fii)* ko + Sh{k,-fii 



ko -£p{k,fiiy 
Zh{k, -fit 



+ F{k*,k,fiiy - F{ko,k,fii) 



and 



p_(/co,k,/ii) = [p+{-k°*,k,-iii)]* . 



(28) 



(29) 



(30) 



Here, Sp{k,fii) and —Sh{k, —fii)* are the two (complex) roots (in ko) of the equation: 



= ko - k + -D+(/co, k, fii) + it 



(31) 
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where iD±{ko, k, fii) are contributions to the inverse, retarded propagator proportional to 
(70 T7'k)/2 arising from interactions. The function F{kQ, k, fii) gives the non-pole part 
of the propagator, and k = |k|. We find that the resonant contributions to the particle 
number-changing sources arise from the pole parts of the propagators, so from here on we 
neglect the terms containing F{kQ, k,fii). 

In the limit g 0, one has Zh and Zp 1, recovering the form of the propagator 
given in Eq. ()25p. For nonzero g, however, Zh is not of order g^ since the particle and 
hole modes arise from mixtures of degenerate states. In particular, at A; = one has 
Zp = Zfi = 1/2. As k becomes large (of order the thermal mass or larger), Z^/Zp <ti 1, 
and the particle dispersion relation is well-approximated by Sp = \kf + m'^{T,fii), where 
m(T, Hi) is the thermal mass. In our particular application to the MSSM, the gaugino Mj 
masses will typically be taken to be of order several hundred GeV, and for the SU(2)2,xU(l)y 
sector, thermal effects do not induce substantial mass corrections. We find that the gaugino 
contributions to the RHS of Eqs. (jl71 118|) are dominated by momenta of order Mj, so that 
the hole contributions to the gaugino S^{x,y) can be neglected. In contrast, for quarks 
we find non-negligible contributions from the low-momentum region, so we retain the full 
structure given by Eqs. ()28H30p in computing their contributions. 

It has been noted in previous studies of quark dampingrates that the one-loop thermal 
widths Tp^h = ^^£p,hik, fi) are gauge-dependent (see Ref. [la] and Ref. [3] therein), whereas 
the thermal masses m(T,fi) entering Sp^h are gauge-independent to this order. Gauge- 
independent widths can be obtained by performing an appropriate resummation of hard 
thermal loops (HTLs) Hi The latter are associated with momenta fco, /c ~ gT, for 

which the one-loop functions D±{ko, k, fi) are of the same order in g as the tree-level inverse 
propagators. In what follows, we will estimate the widths Tp^h based on existing compu- 
tations of damping j2ll . li^ 23l |. deferring a complete computation of the gauge- invariant. 



/ij-dependent contributions in the MSSM to a future study. In general, the residues Zp^h also 
carry a gauge- dependence, and at this time we are not aware of any HTL resummation that 
could eliminate this dependence. In principle, elimination of this gauge-dependence requires 
inclusion of one-loop vertex corrections in the computation of the S'*'(x, ?/) and S''^(x, ?/) ap- 
pearing on the RHS of Eqs. (fT7| ITH|) . and we again defer a complete one- loop computation 
to a future study. 

The derivation of the finite-density scalar propagators proceeds along similar lines. Start- 
ing from the mode expansion of the free scalar field (j){x) in terms of plane-wave solutions to 
the Klein-Gordon equation and following analogous arguments as for fermions, one arrives 
at the following scalar Green's functions: 

G\x, y) = I ^e-^'= (--^)^^(fco, /^.)p(fco, k) (32) 

where the equilibrium distribution functions are: 

gB{uj,fi) = l + nsiuj - fii) (33a) 
93(^,1^) = nsiuj - fii) , (33b) 

with nB{x) = l/(e^/^ — 1) and p(/co,k) given by Eq. ^I^ . As with fermions, one may include 
the effect of thermal masses and widths by replacing — ^ m?{T, /ij) and ie — >■ ie + iT{T, /ij). 
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III. SOURCE TERMS FOR QUANTUM TRANSPORT 



The expressions for G^{x, y) and 5''*' (a;, y) now allow us to compute the perturbative con- 
tributions to the source terms on the RHS of Eqs. ()17I18|) starting from a given electroweak 
model Lagrangian. Here, we work within the MSSM as an illustrative case, but emphasize 
that the methods are general. The Feynman rules giving the relevant interaction vertices in 
the MSSM are taken from Ref. ^|, and in what follows, we only write down those relevant 



for the computations undertaken here. It is useful, however, to place our calculation in a 
broader context by considering the various classes of graphs that generate different terms 
in the QTEs. The simplest topologies are those involving scattering of particles and their 
superpartners from the spacetime varying Higgs vevs (generically denoted v) in the plasma 
[Fig. [T]. These graphs give rise to both the CP- violating source terms discussed in Ref. ^ 
as well as terms proportional to chiral charge. The latter involve the number densities of at 
most two different species, such as the left- and right-handed top quarks [Fig. d^a)] or their 
superpartners [Fig. ^b)]. For purposes of illustration, we follow Ref. (8| and work in a basis 
of mass eigenstates in the unbroken phase, treating the interactions with the Higgs vevs per- 
turbatively. This approximation should be reasonable near the phase transition boundary, 
where both the vevs and their rate of change are small, but it clearly breaks down farther 
inside the bubble wall, where the vevs become large (of order the phase transition temper- 
ature, Tc). In general, one would like to perform a resummation to all orders in the vevs, 
possibly employing the approximation scheme proposed in Refs. 0, 2^. We postpone a 
treatment of this resummation to a future study^. 

Yukawa interactions involving quarks (squarks) and Higgs (Higgsinos) are illustrated in 
Fig. 121 (the self-energies T,^{x,y) are obtained by amputating the external legs). These 



v{x) 

® 



(X) 



v{x) 



® 




(a) (&) 

FIG. 1: Contributions to the relevant self-energies from scattering of particles from the spacetime 
varying Higgs vevs. 



^ The authors of Ref. 11] find that carrying out such a resummation reduces the resonant enhancements 
of the CP-violating sources, but they did not consider the CP-conserving, chirality-changing terms that 
are our focus here. The consistency of the proposed approximate resummation with our power counting 
remains to be analyzed. 



12 



H 



H 





FIG. 2: Contributions to the relevant self-energies from Yukawa interactions 

interactions cause transitions such as f ^ fH, f fH, and / fH. Contributions 
from gauge interactions appear in Fig. El The latter induce transitions of the type / fV, 
f ^ fV, and / ^ fV. In general, one expects the Yukawa and gauge interactions involving 
three different species to depend on sums and differences of the corresponding chemical 
potentials, as in fij — fij — fiy ior the supergauge interactions. In previous studies, it has been 

assumed that the gauginos V are sufficiently light and the coefficients of the corresponding 
terms in the QTEs sufficiently large than one has /ly ~ and /i/ ~ fi^ Although the 
quantitative validity of this assumption could be explored using our framework here, we 
defer that analysis to a future study and take /Xy ~ 0, /x/ ~ /xj. Consequently, one may, as 
in Ref. define a common chemical potential for SM particles (including the two Higgs 
doublets) and their superpartners. 

In previous studies, it has also been assumed - based largely on simple estimates (see, 
e.g., Ref. - that the Yukawa interactions of Fig.|21are sufficiently fast that they decouple 
from the set of QTEs, leading to relations between the chemical potentials for the Higgs 
(Higgsino) fields and those for matter fields. For example, Yukawa interactions that couple 
the Higgs doublet fields H with those of the third generation SU(2)i doublet quarks, Q with 
the singlet top quark supermultiplet field, T, generate terms of the form: 

Fy {hq - Ht + /in) • (34) 

To the extent that Fy is much larger than the other transport coefficients appearing in Eqs. 
(I17I18|) . one has jiq = jiT — f^H plus terms of 0(1/Fy). The remaining terms in the QTEs 
will involve the CP- violating sources, sphaleron terms, and terms that couple left- and right- 
handed chiral charges, such as FM(/iQ — /^t)- Again, this assumption could be tested using 
the current framework, but the computation of Fy is considerably more arduous than those 
discussed below, where we focus on the CP- violating sources and the F^-type terms that 
are generated by the diagrams in Fig. ^ 
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f^H f^H ' f,H f,H 

FIG. 3: Representative contributions to self-energies from (super)gauge interactions. 
A. Bosons 

We consider first the scalar interactions in Fig. ^a). The largest contributions involve 
the L and R top squarks, tL,R owing to their large Yukawa coupling, yt. In the basis of weak 
eigenstates, the relevant interaction Lagrangian is: 

^ = yttLt*R{AtVu - fJ'*Vd) + h.c. , (35) 

where Vu,d are the vevs of -ff^^, and we take v = a/w ^ + f ^ and tan (3 = Vu/vd- Note that in 
Eq. we allow the Vu,d to be spacetime-dependent. In the region of broken electroweak 
symmetry and stable vevs, we have rrit = ytVu- 

Using the Feynman rules for path ordered perturbation theory, it is straightforward to 
show that the diagrams in Fig. ^a) generate contributions to T,ji{x,y) of the form: 

^r{x, y) = -g{x, y)Gl{x, y) (36) 

where 

9{.x, y) = yf [AtVuix) - iJ,*Vdix)] [Alvu{y) - i^Vdiy)] . (37) 
Substituting Eq. ^ into Eq. jTTj) leads to: 

d.i^ni^) = S,^{x) (38) 

for right-handed top squarks, where is the corresponding current density and the source 
Si^{x) is 



j (fz j dzo !^[g{x, z) + g{z,x)]Re[Gl{x, z)G'^{z,x) - Gt{x, z)G^{z,x)] 
+i[g{x, z) — g{z, x)] Im[G'^(x, z)G'^{z, x) — Gf{x, z)G'^{z, x 



(39) 
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where the L,R subscripts indicate the propagators for the L and R top squarks. 

The first term in the integrand of Sf:^{x) is CP-conserving and leads to the FM-type terms 
discussed above, while the second term in the integrand provides the CP- violating sources. 
We concentrate first on the former. Expanding g{x, z) about 2; = x it is straightforward to 
show that only terms involving even powers of derivatives survive in g{x^ z) + g{z, x). Under 
the assumptions of gentle spacetime dependence of the Vi{x) near the phase boundary, we 
will neglect terms beyond leading order and take g{x, z) + g{z^ x) ~ Igi^x, x). Consequently, 
the CP-conserving source is: 

d^z I dzQ [Gl{x,z)G^{z,x) -Gt{x,z)G'^{z,x)] (40) 
J —00 



-00 

dh j dzoj j^^j ^e-('=-'^)-(--)pi(A:o,k)p«(go,q) 
X [9B{ko,^^L)9^{qo,^^R) - ^5(^0, /^L)^s(go, /^i?)] , 



with 



g{x, x) = [\fi\hj{x) + \At\'^vl{x) - 2vd{x)vu{x) Re(pAi)] . (41) 

Performing the d^z integral leads to a 5 function in momentum space. After carrying out 
the d^q integral, we perform the contour integrals for feg and Qq, expand to first order in 
fiL^n/T, and obtain: 



cCP/ N -'- ^cUt I A /■ \ * / N|2 f°° dkk"^ 



~ £l + £r 



(42) 



where 



a;i,^=|k|^ + M^^^ (43a) 

^L,R = ^^L,R - i'rL,R (43b) 

= - (ex/T_i)2 ' (43c) 

and Mf^ ^, Tl r are the thermal masses and widths for the ti^ij, and the factor of Nc comes 
from summing over the colors. Note that, in arriving at Eq. (02)), we have neglected 
the /ij-dependence of the pole residues Z(T, fiL^n), thermal frequencies, u!l,r(T, fiL^u), and 
widths, TL^Fi(T, iil,r)- The effect on S^^{x) of the /Xj-dependence of the residues and thermal 
frequencies is sub-leading in the gauge and Yukawa couplings, whereas the effect from the 
thermal widths occurs at leading order. The /ij-dependence of T l^r{T, iil,r) is simply not 
known, however, so we do not include it here. A more explicit expression for the dependence 
of Sf^ix) on the thermal frequencies and widths is given in Eqs. (IB1IIB3I) Appendix iBl 
For purposes of future analysis, it is useful to rewrite Eq. as: 



where 



Sf^ = Tf (/i^ + /xh) + Fr (p^ - p^) , (44) 



t y 4^2 I ^ "V / A- «v n ^^^^ ^ Sl-S*^ Sl + Sr 

(45) 
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Before proceeding with the CP- violating source, we comment briefly on the structure of 
Eqs. ()44ll45p . In particular, we note that 

(i) Terms of the type are absent from the conventional QTEs for EWB. It is straight- 
forward to see that in the absence of interactions that distinguish between II and t^, 

= 0, as the integrand of Eq. ()45p is antisymmetric under L ^ R interchange. In 
contrast, the transport coefficient is nonzero in the limit of exact ^ tji symme- 
try. This term corresponds to the usual damping term in the QTEs associated with 
scattering from the Higgs vevs. 

(ii) In the absence of thermal widths r^ ij, the quantity in brackets in Eq. is purely 
real, and so the damping term would be zero. 

(iii) The structure of the energy denominators implies a resonant enhancement of the 
integrand for ~ . A similar effect was observed to occur for the CP- violating 



sources (see below) in Refs. j8|, |26|. The expression in Eq. makes it clear that 
the relaxation terms display a resonant behavior as well. The resulting quantitative 
impact of this resonance baryon asymmetry is discussed in Sect. IVl 

Properties (ii) and (iii) are shared by all source and damping terms, we discuss below. Note 
that the explicit factors of /iL,_R/T and property (ii) imply that, away from the resonance 
region, S^^ is C(e^). 

The computation of the CP- violating source, given by the second term in Eq. ()39p. 
proceeds along similar lines. In this case, the coefficient [g{x, z) — g{z, x)] vanishes for x = z, 
so we must retain terms at least to first order in the expansion about x = z: 

g{x,z) - g{z,x) = 2iy^ lm{ij,At) [vdix)vuiz) - Vd{z)v^{x)] 

= 2iy^ lm{^At){z - x) [vd{x)dxVu{x) - Vu{x)dxVd{x)] H , 

where the indicate higher order terms in the derivative expansion that we neglect 

for the same reasons as discussed previously. When the linear term in Eq. (|46p is sub- 
stituted in Eq. (|39p . only the time component yields a nonzero contribution. The spatial 
components vanish due to the spatial isotropy of the spectral density: g^^ko, fi)p{ko,k) = 
g'^{ko, fi)p{kQ, |k|). We may then make the replacement: 

g{x, z) - g{z, x) 2iy^ lm{fiAt) [vd{x)vu{x) - Vd{x)vu{x)] {z - xf 

(47) 

= 2iyl \m{pAt)v{xf(3{x) {z - xf . 

In general, we expect /5 to be of order Vwlc, so that the CP- violating source is first-order in 
one of the small expansion parameters discussed earlier. Consequently, when evaluating this 
term, we may neglect the /iL^^-dependence of the g^{ko,fi). After carrying out the {ko,qo) 
contour integrals and performing the time integration, we obtain: 

S^-^lra(uAAv(x?3(x) H dkk^ J usjS^n) - nsiSL) I + n.jSn) + nsjS,) ] 

(48) 

Again, property (ii), in conjunction with the factor of /5 oc Vw, implies that S*?^ is 0{e^). 
An expression giving a more explicit dependence on the widths and frequencies appears 
Eq. ()B4|) of Appendix ^ We note that our result agrees with that of Ref. except for a 
different relative sign in front of the cos 20 term of that equation and the overall factor of 
Nc. 
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B. Massive fermions 



The computations for fermions proceed along similar lines. We consider first the source 
terms for Higgsinos. We recall that it is useful to redefine the Higgsino fields to remove the 
complex phase from the Higgsino mass term: 



/rnass 



via 



leading to: 



Defining the four component spinors, 



for the Higgsinos, and 




Hi 



HO 



-HI 
H'J 



(49) 
(50) 
(51) 

(52) 



wo 




^ B 



B 

5t 



(53) 



for the gauginos, leads to the Higgsino-gaugino-vev interaction: 



1 , 



(54) 



Note that the spinors ^'^o and jj+ satisfy a Dirac equation with Dirac mass even 
though the H^^ are Majorana particles. The '^^± are Dirac particles of mass M2, whereas 
the \Ef^o and ^'^o are Majorana particles with Majorana masses M2 and Mi, respectively. 
We also note that the construction of the Dirac spinor ^ fjo allows one to define a vector 
charge and corresponding chemical potential, /i^o, for the neutral Higgsinos, even though 
they are Majorana particles. In contrast, there exists no such vector charge for the '^'^0 and 
ipfjo- One may, however, study the quantum transport of the axial charge of the Majorana 
fermions using Eq. ()2()|1 . An attempt to do so for the neutral Higgsinos was made in Ref. 



IjJ, though only the CP- violating sources were evaluated using non-equilibrium methods. 



The impact of the corresponding axial charge density on the baryon asymmetry was found 
to be small. We return to this issue in a future study, and consider only the vector densities 
below. _ 

The most straightforward computation is that of the if^ source terms. For notational 
convenience, we rewrite the chargino interactions in Eq. ()54|1 as: 



-92^ H+ [9l{x)Pl + gR{x)PR] + h.c. 
In this case, the self-energy generated by Fig. HJ^a) is: 

T,^±{x,y) = -gl [gL{x)PL + 9r{x)Pr] Sy^±{x,y) [gL{y)*PR + 9R{y)*PL 



(55) 



(56) 
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Defining: 



9A{x,y) = Y [gMgLivY + gRix)gR{yy] 
gB{x,y) = Y [gL{x)gRiyy + gRix)gLiyy] 



(57a) 
(57b) 



we obtain for tlie RHS of Eq. (ITH|l : 



! Sz f dzo J2 I 

J J -co j=A,B ^ 



(58) 



[gj{x,z)+ gj{z,x)] Re Tr S^^{x, z)S^^{z, x) - S^^{x, z)S^^{z, x) 



+ i[gj{x,z) - gj{z,x)] Im Tr S^^{x, z)S^^{z,x) - S^^^{x, z)S<;^^{z,x) 



where the subscripts "A" and "B" on the traces denote the contributions arising from the 
}i and m terms, respectively, in the spectral function in Eq. (an overall factor of 1/2 
due to the presence of the chiral projectors Pl,r has been absorbed in the definition of the 
gA,B)- 

As in the case of the scalar fields, the leading density- dependent, CP-conserving contri- 
bution to S^± (x) arises from the term in Eq. containing the x ^ z symmetric factors 
[gj{x,z) + gj{z,x)]. To lowest order in Vw, we may set x = 2; in these factors. Using the 
spectral representation of the S^{x,y) given in Eq. ()25p. including gauge- invariant thermal 
masses and widths, and expanding to first order in /ij/T, we obtain the chirality-changing 
source term: 

where 



H' 



Im<^ M2 cos6'^sm2/5J 



^ii^w 



[Sy^S^ + k'^ - M2 Ifil cos 6'^ sin 2(3] 



H) 



(60) 



where the definitions of o;^^ and £jj^ are analogous to those given in Eqs. and 

hpi^x) 



(61) 



Also, the factor of cos 6^ is very nearly 1 for the region of small 6^ in which we find ourselves 
in subsequent sections. The explicit dependence of on thermal frequencies and widths 
is given in Eq. ()B5|1 of Appendix iBl 

In the present case, we follow Ref. and assume no net density of gauginos, thereby 
setting /i^^± = in Eq. (jSHI) and giving: 



(62) 
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with r^± = rt^ + r^^. in this case, it is straightforward to obtain the corresponding 
source term for the neutral Higgsinos, 

where F can be obtained from the formulae for F ^± by making the following replacements: 
92 -^92/ \/2 for intermediate states and (72 9i/ v^, i^^y ~^ ^^'i Tjj for the 

B intermediate states. 

The Higgsino CP- violating source arises from the second term in Eq. As before, 

we expand the gj (x, z) to first order about x = z and observe that only the Xq — yo compo- 
nent survives when the d^z integration is performed. Also note that gA{x,y) — g^iVyX) = 
2i lmgj^{x,y) = so that only the terms proportional to the Higgsino and gaugino masses 
contribute. The result is: 

sfA^) = ^v{xf^{x)M, |/i| sin 5^ 

r-^dk^ / nH%)-n^(£:t) 1 - - nHg^) 1 ^^^^ 

Jo c^hc^^ "^i (%-^~)^ i^w+^ny )' 

The corresponding expression for S^{x) can be obtained by making the same replacements 
as inclicated above for the CP-conserving terms. The correspondence with the results of 
Ref. 9] can be seen from Eq. ()B7|) of Appendix^ We again find essential agreement, apart 
from a sign difference on the cos 20 term. 



C. Chiral fermions 

The final source term associated with Fig. C^a) involves L and R top quarks. At this 
order, the latter only contribute a /Xj-dependent CP-conserving term. In order to illustrate 
the structure of this term that arises when the terms of 0{g'^) are retained, we employ the 
interacting fermion propagators of Eqs. ()28II30|) . The result is: 

S^ai^) = (/^iL +f^tn)+ Tf" ifitL - /iifl) , (65) 

with 

Z^(k)Z^(k)\ , „ , 1 

+ sL _ s J [h^i^i) T h,{St)] + (p - /^) I . 

Here, the "p" and "/i" subscripts indicate contributions from the particle and hole modes, 
and "L " and "P" refer to left- and right-handed quarks. We have not included in our 
calculation the effects of /Xi^ ^-dependence of the widths Fp'^(T, /Xf^ ^), which in principle 
also enter at this order. For an expanded version of Eq. (jHSl), including these effects, see Eq. 
finsl) in Appendix IB 
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In the limit of ^ ti? symmetry, Tf^ vanishes, and T^^ simphfies to: 

1 NcvXixf r,,,2. f ZM\ , ZH{k)\ . 

= 7f 2 / dkk lm< — — hF{Sp) + — — hF{Sh) 

J- Jo K . „s 

We observe that all contributions to the CP- violating source terms and the vanish in the 
limit of zero thermal widths. Since the widths are generically of order g'^T (here, g denotes 
either a gauge or Yukawa coupling), the source terms for the QTEs are generally fourth 
order in the couplings. 



IV. QUANTUM TRANSPORT EQUATIONS AND ps 

We now discuss diffusion equations for the particle species that significantly contribute 
to the density of left-handed doublet fermions ni^x) [cf. Eq. (j2))] that acts as the "seed" for 
baryogenesis. We subsequently relate pb to and solve explicitly the equations in the case 
of a simple geometry and profile for the bubble wall describing the phase boundary. 



A. Solving the diffusion equations 

Using the source terms computed in Section 11111 one can arrive at a coupled set of dif- 
ferential equations for the various particle number densities. These equations simplify con- 
siderably under the assumptions of approximate chemical equilibrium between SM particles 
and their superpartners (/i/ ~ fij with /x^ pa 0), as well as the between different members of 
left-handed fermion doublets {pw^ ~ 0). In this case, one obtains transport equations for 
densities associated with different members of a supermultiplet. This approach is the one 
followed in Ref. [13], and for pedagogical purposes we summarize the development here. 

First, we define the appropriate supermultiplet densities: 

Q = nt^+ ni^ + rib^ + (68a) 

T = nt^+ ni^ (68b) 

5 = rife^ + n^^ (68c) 

H = njj++ Uho - riH- - riHO + nfj+ - n^- + rifj^ - n^o , (68d) 

where the Higgsino densities arise from the vector charges = ^^^+7°\l/^+ and n^o = 
^fj^o7°^/^o associated with the Dirac fields defined in Eq. (|^. There are analogous defini- 
tions for the first- and second-generation (s)quarks. Although we do not consider them here, 
one may also define the corresponding axial charge densities. In the case of the Higgsinos, 
for example, it will involve the sum, rather than the difference, of the u- and (i-type Higgsino 
densities^ 

^ This density was considered in Ref. ^'^'^ its overall impact on the baryon asymmetry found to be 

small. 
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The diffusion equation for a density rii has the structure: 



d,Jt = Sr + Sf + (69) 

where Jf is the current associated with the density rii, Sf^ and Sf^ are the source terms 
computed above, and Sl"^^ is the sphaleron transition term. Various derivations of the strong 
sphaleron term appear in the hterature, so we do not reproduce them here. However, we 
note that the expressions in Refs. [H |27| have erroneously omitted a factor of 1/Nc '28']. 

The CP-conserving damping terms Sf^ have been given in Eqs. (HH), (j^ . and (jU^ to 
hnear order in the appropriate chemical potentials. Assuming local thermal equilibrium we 
relate the number densities to the chemical potentials via: 

ni = J [N{ujk, Hi) - N{ujk, -fii)] , (70) 

where N{uj,h) is the appropriate boson or fermion distribution function and Qi counts the 
internal degrees of freedom (spin and color). Dropping terms of 0{fif), one obtains: 

rii = /ii , (71) 

o 

where the factors ki{mi/T) are exponentially small in the regime m^/T ^ 1, and reduce 
in the massless limit to /^^(O) = 1 for chiral fermions, /cj(0) = 2 for Dirac fermions, and 
/ci(0) = 2 for complex scalars. In our analysis we keep the full dependence on irii/T: 

ki{m,/T) = h{0)^ / dxx . ,^,,, V^'-myT\ (72) 

^ Jm/T ± 

where for fermions (bosons) cf(^b) = 6 (3), and we choose the +(— ) sign in the denominator. 
Using Eq. ((ZH) in Eqs. ()44l59iei5jl . and defining: 

= ^ (r,^ + rf) (73a) 
r. = ^ (r^± + r^o) , (73b) 

the resulting set of coupled transport equations is: 



kr kn kn J V kn kx ks 



t 



a.<3,.-ri(£ + ^)+r„(£-|) (T4b) 



ykr kn kq J ykq kx k 



t 



^"^^-r.-.r^^^-.-USr. (74c) 
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where Tgs = Gn'^ajT, with k' ~ 0(1). 

In writing down Eqs. ()74all74c|) . we have also included the Hq^qji and Hq^q^ Yukawa 
interaction term that arises from Fig. 2, though we have not computed the corresponding 
transport coefficient Fy. In this work we will again follow the authors of Ref. |lO|, who 
estimate Fy ^ F^j. For k' ~ 0{1), one also has F^^ ^ F^. These facts allow one to alge- 
braically relate the densities Q and T to if, by setting the linear combinations multiplying 
Fy and Vgs equal to 5y = 0(1/Fy) and 5ss = 0{1/Vss), respectively. One then obtains: 

{kB - 9kT)kQ 7T , r , r fy^ \ 

^ = i9kr + 9kQ + ks)k^ ^ + + 

{9kT + 2kB)kT ^, ^ ^ ^ 

with known coefficients aQY,Qs,TY,Ts- Taking 2 x [Eq. (ITiall ] + [Eq. (ITiEll ] + [Eq. (I7ic|l ]. 
introducing the diffusion approximations T = —Dg'VT, Q = —Dg'VQ, H = —Dh'VH, and 
using Eq. (f7H|l leads to: 



H - DV'^H + TH - S = 0{5ss,Sy) , (76) 



where^ 



^ ^ (9/cQfcr + fc^feg + 4:kTkB)Dq + kH{9kT + 9kQ + kB)D 

9kQkT + kBkg + AkrkB + kH{9kT + 9kQ + k ' 
- ^ {9kQ + 9kT + kB){T-M + Th) - {SkB + 9kQ - 9kj 
9kQkT + kBkq + AkrkB + kH{9kT + 9kQ + k 



■b) 



(77a) 
(77b) 



kH{9kQ + 9kT + kB) ( qQ^^ ^ qQ^^ 

9kQkT + kBkq + AkrkB + A;j^(9/ct + 9kQ + A;b) V 



The subleading terms 5y,ss can be determined by use of Eqs. (ffH|) in Eqs. ()74a|74b|) . We 
include the effect of 5ss in our final expression for pB [10], although its effect is negligible in 
the relevant MSSM parameter region. 

Equation (fTUI) can now be solved for a given set of assumptions about the geometry of 
the bubble wall. Again, for clarity of illustration, we will work in a framework that allows 
us to carry analytic calculations as far as possible, leaving to the future a numerical solution 
of the equations for a realistic wall geometry and profile. First, as commonly done in earlier 
studies, we ignore the wall curvature, thereby reducing the problem to a one-dimensional 
one in which all relevant functions depend on the variable ^ = |x + v^t|, where is the 
wall velocity. Thus, 2 < is associated with the unbroken phase, ^ > with the broken 
phase, and the boundary wall extends over < ^ < L^. Second, we take the relaxation 
term F to be nonzero and constant for ^ > 0. The resulting solution for H in the unbroken 
phase z < (related to pb as shown below) is: 

H{z) = Ae''-~''^ (78) 



^ Our expressions differ from those in Ref. IQ], which we beheve result from an algebraic error. The 
numerical impact of this difference, however, is not significant. 
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with 



The above equation is vahd for any shape of the source S{z). For simphcity, however, we 
assume a simple step-function type behavior for the source: S nonzero and constant for 
< z < Lw Speciahzing to this case of constant sources in < z < L^, using 4Z)f ^ v^^, 
Lw^/tJd ^ 1, and taking f = rr (Th + r^_^) from Eq. ()77a|) . we arrive at: 

A = kHL^\l^ f *_ = • (80) 




When evaluating the source terms S'^, S?^ [see Eqs. fl^ . ljMj) ] for this simple profile one has 

to use (3 = Vj^Ap/L^: thus A is explicitly proportional to and is only weakly dependent 
on Lyj. Solutions for Q and T are then obtained via Eqs. (f75j) and (f?^. 



B. The baryon density ps 



Neglecting the wall curvature and assuming a step- function profile for the weak sphaleron 
rate, the baryon density satisfies the equation jl^. l29l|: 

DMz) - vM^) - 0{-z)npB = e{-z) ^r„,ni(z) , (8i) 

where ni? is the number of fermion families and the relaxation term is given by j2^: 

^ = rws[?fi + ^)" + ^l , (82) 



where rigq indicates the number of light squark flavors, and the weak sphaleron rate is given 
by = 6fi;a^T, with k ~ 20 (20]- 

The solution to Eq. (|HT| in the broken phase, eventually growing into the Universe, is 
constant and given by: 



PB 



2f,„ 



nL{x) e"^/"- dx . 



(83) 



Neglecting leptonic contributions, n^is given in the unbroken phase by the sum of left-handed 
quark densities over the three generations (Qil, Q)- Since appreciable densities of first 
and second generation quarks are only generated via strong sphaleron processes, it is possible 
to express Qil and Q2L in terms of Q and T, in such a way that ul = Q + Qn + Q2L = 
5Q + 4T ilO|. Using then Eq. ^ one obtains : 



-H 



ri + r2 



D 



(84) 
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where 



ri 



r2 



kH{9kQ + 9kT + ks) 
k%{5kQ + AkT){kQ + 2kT) 
kH{9kQ + 9kT + kBY 



(85a) 
(85b) 



and finally, in the broken phase: 
Pb{z > 0) 
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(86) 



where the second line is true in the limit ^ ADgTZ, which holds for the parameters we 
have chosen in this calculation. The contribution from the first term in Eq. (j86p is linear in 
Viu, due to the linear dependence on contained in the /3 appearing in the CP- violating 
sources. The second term is suppressed by two additional powers of and generally leads to 
a negligible contribution to in the MSSM case (see discussion below). It could, however, 
be dominant in the case of heavy degenerate ii and tj^, which leads to ri ~ |lO|. 

The central feature emerging from the above discussion is that the net baryon density 
is proportional to ^ ~ S^/VT. A large relaxation rate T for the relevant charges will 
suppress the overall baryon asymmetry. While in Refs. |2^ it was pointed out how a non- 
equilibrium quantum transport could result in a resonant enhancement of we observe 
here that similar resonance effects in the relaxation terms will mitigate the impact of the 
enhanced sources. In the next section we discuss the numerical impact within the MSSM, 
but caution that reaching definitive conclusions will require computing the other transport 
coefficients, such as Fy, within the same framework. 



V. BARYOGENESIS AND ELECTROWEAK PHENOMENOLOGY WITHIN 
THE MSSM 

The results derived in the previous Sections allow us to perform an illustrative, prelim- 
inary analysis of baryogenesis within the MSSM. This should be taken as an exploration, 
whose robustness will be tested once we implement the next steps in our treatment of the 
source terms in the transport equations. With this caveat in mind, we explore the connec- 
tions between electroweak baryogenesis and phenomenology within the MSSM, focusing in 
particular on the implications for EDM searches. Throughout, we assume that all the terms 
in the Higgs scalar potential and all gaugino masses are real, while all the A-parameters 
(trilinear scalar couplings) are equal at the GUT scale, therefore sharing the same phase 
(pA- In this case, the baryon asymmetry and EDMs are sensitive to the two independent CP 
violating phases 0^ and 0^- 

From the structure of Eqs. ()86l80p and ()48I64|) we can write the baryon-to-entropy density 
ratio"^ Yb = Pb/s as: 

Yb = Fi sin0^ + F2 sin (0^ + 0a) , (87) 



We evaluate the entropy density at the electroweak phase transitions via s — {2n^)/A5 x geff{T) T^, with 
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FIG. 4: Left panel: CP-violating Higgsino source Sj^ = /(v^/Ssinc/)^), as a function of 

Right panel: relaxation rate Ry = {Th + T~^^) / iVh + '^m)H.N.-, normalized to the value used in |lo| . 
as a function of We have taken M2 = 200 GeV, and the values of all other parameters as 
indicated in the text. 



where we have isolated the dependence on the phases 0^ and 0^. The first term that contains 
Fi stems from the Higgsino source, while the F2 term arises from the squark source. 

The functions Fi and F2 display a common overall dependence on bubble wall parameters 
(v^, Lyj, A/9), while having distinct dependence on other MSSM mass paramters such as 
the soft mass parameters for gauginos (Mi^2) and squarks (Mf^, M^^), the triscalar coupling 
\At\, and tan/3. In order to assess the size of Yb and the impact on CP- violating phases, we 
must choose a reference region in the MSSM parameter space, and we follow two obvious 
guidelines: (i) we require that v{Tc)/Tc ^ 1, so that the baryon asymmetry is not washed 
out in the broken symmetry phase; (ii) we require no conflict with precision electroweak 
physics and direct collider searches. Both criteria lead to non-trivial restrictions. 

The condition of a strongly first-order phase transition [v{Tc)/Tc ^1] requires light scalar 
degrees of freedom coupling to the Higgs sector. It has been shown [3l|, |3j| that within the 
MSSM the only viable candidate is a light top quark, which should be mainly right-handed 
(tR) in order to avoid large contributions to the p parameter. Quantitatively, for lightest 
Higgs boson mass rrih ^ 120 GeV, one needs 100 GeV ^ < mt, and sufficiently small 
stop mixing parameter \At — ji/ tan/5| < 0.6Mf^ [sH. Moreover, present experimental limits 
on rrih and the constraint v(Tc)/Tc > 1 jointly require either values of tan/3 > 5 or M^^ in 
the multi-TeV region Based on these considerations, for illustrative purposes we work 
with the following values of MSSM parameters at the electroweak scale: M^^ = 0, M^^ = 1 
TeV, \At\= 200 GeV, M2 = 200 GeV, tan /? = 10. We also take for the CP-odd Higgs mass 
rriA = 150 GeV, which translates into A/? ~ 0.015 0]. We vary in the plots the scale in 
order to display the resonant behavior for ~ M2. Finally, for the bubble wall parameters 
we adopt the central values = 0.05 and = 25/T [33[. 

With the above choice of MSSM parameters, the stop-induced contribution to Yb is 
suppressed (P2 ~ lO^^Pi), since one is far off the squark resonance [(M^^ — Mf^) ^ ^t^)]- 



gcs — 130.75, resulting in s ~ 57.35 T^. Similarly, to convert the present ratio pB/n-y to Yb, we use the 
relation s = 7.04 n^. 
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FIG. 5: Left panel: Higgsino contribution to Yb { CfE,q. ()87|) '). normalized to the observed value. Fi 
displays residual resonant behavior for ~ M2. All other input parameters are given in the text. 
Right panel: Stop contribution to Yb (C/Eq. (|S7j) ) normalized to the observed value. The upper 
curve is for Mr = Mr , while the lower one is for Mr 3> Mr . We have taken here Mr = 100 
GeV, \fi\ = 200 GeV, and have allowed M^^ to reach unrealistically low values to explore the size 
of the squark resonance. For realistic input parameters F2 <C -Fi. 

On the other hand, the Higgsino-induced contribution Fi can account for the observed Yb 
even without maximal values of | sin0^|. We highlight below the salient results of our study: 

• The primary result of our analysis is that both the source and the relaxation term 
Th display the resonant behavior jl,|2^ typical of quantum transport for ~ M2. We 
illustrate this in Fig. 01 the left panel shows the behavior of the rescaled CP- violating 

higgsino source Sfj = — S'?^/(f sin 0^) versus while the right panel displays the 
ratio Rt of the relaxation term (Th + T^) as calculated in this work to the one used 
in previous studies, (Th + T]^^)h.n. To our knowledge this is the first explicit 

calculation showing resonance behavior for the relaxation term F ~ r^iTh + T^). 

• Since Fi is proportional to S'^/ ^/I\~+T^, the baryon asymmetry retains a resonant 
behavior, albeit with an attenuation of the peak due to the enhanced relaxation term. 
This is shown explicitly in Fig. El In the left panel we plot Fi/Y^^^^, normalizing to 
the baryon asymmetry extracted from CMB studies [3|: y^^^^p _ 2 ^ 10~^^ 
(the quoted error corresponds to 95% CL). 

• For completeness we also display in Fig. El (right panel) the behavior of the squark 
contribution F2/Y^^^^ as a function of M^^, with M^^ = 100 GeV. Within the 
MSSM, precision electroweak data and the requirement that v{Tc)/Tc 1 force the 
masses to be far away from the peak region. However, in extensions of the MSSM 
where the phase transition is strenghtened by additional scalar degrees of freedom this 
contribution might be important (see, e.g., Refs. j^iS])- 

• For given values of the MSSM parameter space explored here, successful EWB carves 
out a band in | sin(^^| centered at | sin0^| = Fj'^^/lFil (whose width depends on the 
uncertainty in Y^^^). Due to the resonant behavior of Fi, the location of this band 
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FIG. 6: Allowed band in the | sin (/)^|-|//| plane, obtained by requiring successful electroweak 
baryogenesis. All other MSSM parameters are given in the text. The light-shaded (green) narrow 
band corresponds to the experimental input from WMAP, while the two bands combined [dark 
(blue) + light (green)] correspond to input from Big Bang Nucleosynthesis. 



is highly sensitive to the relative size of M2 and As illustration, in Fig. IHl we 
plot the allowed band in the | sin0^|-|/i| plane determined by the baryon asymmetry, 
with all other MSSM parameters fixed as above. The bands in the plot combined 
together correspond to the baryon density determined from Big Bang Nucleosynthesis, 
yBBN ^ (7.3±2.5) X 10"" (the error corresponds to 95% CL 0). Using WMAP input 
leads to the narrow, lighter-shaded band in our plot located at the upper edge of the 
BBN-induced band. 

We conclude this section with a brief account of the connections between the baryon 
asymmetry and EDM phenomenology. Since the Standard Model predictions are in general 
highly suppressed and well below present experimental sensitivity, limits on the electron, 
neutron, and atomic EDMs can be used to constrain the phases of a given new physics 
model. Present limits of interest to us are: 

27„ /nKCW nj \ Jo2 I ^ o 1 X. 1 n-28 



141 < 1.9 X 10~^'e ■ cm (95% CL) j3| , \dHg\ < 2.1 x 10~^«e ■ cm (95% CL) [3 

Although a single EDM can be sufficiently small even for maximally large CP-violating 
phases (du e to cancellations), constraints from more than one EDM can be very powerful. 
In Ref. m, for example, it was pointed out how limits on electron and ^^^Hg EDMs single 



out a well defined region in the (j)fj,-(f)A plane, for given values of gauginos, squark and 
slepton masses. As shown above, for each point in the MSSM paramter space, electroweak 
baryogenesis also selects a band in the (pfjr'pA plane. This implies in general non-trivial 
constraints on the MSSM parameter space, as the EDM-allowed region need not in general 
coincide with the one required by the baryon asymmetry. 

To illustrate this situation, we have evaluated the bands in 0^ — (pA allowed by present 
limits on electron EDM, mercury EDM, and EWB for two representative points in the MSSM 
parameter space (see Fig. [7j). In our analysis we take the expressions for the electron EDM 
and quark chromo-electric dipole moments from Ref. [s^. In relating the ^^^Hg EDM to 
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FIG. 7: Allowed bands in the (pa^ipA plane implied by consistency with the 



C.L. limits on 



electron EDM (141 < 1-9 x lO'^'^e • cm 0), mercury EDM {\dHg\ < 2.1 x lO'^^e • cm j35|), and 
baryo^nesis. The shaded [dark (blue) and light (green) combined] EWB band corresponds to BBN 
input ;2'], while the narrow light-shaded (green) band on the left corresponds to WMAP input [3]- 
In the left panel we use \fi\ = M2 = 200 GeV (resonance peak), while in the right panel we use 
M2 = 200 GeV and = 250 GeV (off resonance). In both cases the other supersymmetric masses 
are as specified in the text. 



the quark- level CP- violating couplings, we follow the treatment of Ref. |36|] , where it was 
shown that the dominant contribution arises from the chromo-electric dipole moments of 
quarks (dq) according to 



d 



Hg 



dd-du- 0.0124 X 3.2 ■ 10~^e 



(88) 



According to the same authors, the analysis of the neutron EDM involves additional effects, 
such as the CP-violating three-gluon operator CGG, that require a detailed analysis going 
beyond the scope of the present work. Although the experimental bounds on the neutron 
EDM — as well as the prospects for future improvements — are competitive with those for the 
electron and neutral atoms, we defer an analysis of its implications for EWB to a future 
study. We also neglect the renormalization group evolution of 0^ and (pA from the weak scale 
to the atomic scale, having assumed a common, flavor-independent phase for the tri-scalar 
coupling at the former. 

The plots in Fig. [7| correspond to taking the first and second generation sfermions, along 
with the gluinos, to be degenerate with masses equal to 750 GeV; the gaugino mass Mi = 100 
GeV; and the triscalar coupling A = 200 GeV. We consider then two cases for M2 and fi: 



For a recent reanalysis of hadronic EDMs in SUSY see Ref. |3 
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the left panel corresponds to the resonance peak M2 = = 200 GeV, while the right panel 
corresponds to off-resonance parameters M2 = 200 GeV, = 250 GeV. For these choices 
of MSSM parameters, Eq. ()87|) predicts for Yg: 

M2 = = 200 GeV : Yb = -1.3 x 10"* sin + 1.7 x 10"^^ sin((^A + 
M2 = 200 GeV, = 250 GeV : Yb = -2.0 x lO"%in0^ + 4.6 x 10"^^ sin(0A + (p^) 

These cases illustrate the main trend: for M2 ~ |/i| electroweak baryogenesis requires rel- 
atively small phases, and is consistent with the constraints from EDMs. As one moves 
off resonance, then larger phases are needed to generate the observed baryon asymmetry, 
and this requirement tends to conflict with the EDM constraints. Indeed, within our sim- 
plified analysis, we find that baryogenesis and EDM constraints become inconsistent for 

— M2 ^ 50 GeV, when all other superpartners are kept around 750 GeV. Of course, in- 
creasing (decreasing) the sfermion masses relaxes (tightens) the EDM limits on CP-violating 
phases and affects the above conclusion. 

Ultimately, if supersymmetry is discovered at collider experiments, spectroscopy will 
dictate the input for mass parameters. Then joint constraints from low-energy EDM mea- 
surements and collider searches could be used to tightly test the scenario of baryogenesis at 
the electroweak scale. 



VI. CONCLUSIONS 

It is instructive to consider the essential physics leading to the enhanced sources and 
relaxation terms discussed in this work. The propagation of quasiparticles in the plasma 
is modified by scattering from the spacetime varying Higgs vevs that causes transitions to 
intermediate states involving other quasiparticle species. The system retains some memory 
of each scattering due to the presence of thermal widths, Fj, that reflect the degeneracy 
of states in the thermal bath. For Fj = 0, the oscillating exponentials appearing in the 
Green's functions wash out any memory of the scattering. For Fj ^ 0, the Green's functions 
now contain decaying exponentials as well as oscillating terms, and the memory wash out 
is incomplete. The impact of quantum memory effects are, thus, characterized by the ratio 
of time scales, Tint/Tp ~ Ti/ui, where Tint is the characteristic propagation time associated 
with a quasiparticle of frequency Ui and Tp ~ 1/Fj, the plasma time, is time scale on which 
transitions between the quasiparticle and other, degenerate states may occur. To the extent 
that the quasiparticle thermal mass and/or three-momentum is large compared to F,, this 
ratio Tint/Tp is C(e). 

A special situation arises, however, when the spacetime variation of the Higgs vevs is 
gentle and the thermal mass of an intermediate state is close to that of the initial state. 
Under these conditions, the scattering event injects essentially zero four-momentum into the 
initial state i, leading to resonant production of the intermediate state j. The characteristic 
lifetime of the latter is no longer Tint ~ l/co'j, but rather the resonance time scale 

(89) 



Aa;2 + F^^ 



where Au = Ui — ojj and Fj^ = Fj + F^ [see, e.g., Eqs. flBlllB3|) and ()B5IIB7|) of Appendix 
|B]. In this case, the impact of quantum memory is characterized by the ratio Tres/Tp. For 
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|Aa;| <^ Tij, this ratio becomes of 0{1), and the impact of quantum memory is resonantly 
enhanced^. On the other hand, for \Auj\ ^ Tij, the ratio is 0{e) and one returns to the 
more generic conditions. 

In this study, we have shown how this effect can enhance both the particle number- 
changing relaxation terms as well as the CP- violating sources that enter the transport 
equations relevant to electroweak baryogenesis. Importantly, the effect of resonant relaxation 
tends to mitigate the impact of resonantly-enhanced sources, as both enhancements occur 
under the same conditions for the electroweak model parameters (in this case, those of the 
MSSM). We suspect that analogous resonant effects occur in other transport coefficients, 
such as the Ty Yukawa terms discussed above, but that the conditions on model parameters 
leading to enhancements — owing to simple kinematic considerations — will be different. It 
may be, for example, that the Yukawa interactions are no longer fast compared to the 
Higgs vev induced transitions when the latter are resonantly enhanced, and in this case, the 
solution to the differential equations will differ from the general structure obtained here and 
by other authors. This possibility is one that should be explored in future work. 

Additional refinements of the present analysis are clearly in order, including some form of 
all-orders resummation of the Higgs vev insertions (possibly along the lines proposed in Refs. 



ll|,|25|) and a treatment of the axial charge transport equations via Eq. (j^Uj) . In principle, 
one would also like to study the density dependence of the thermal frequencies and widths, 
the impact of nonzero gaugino densities, variations in bubble wall geometry, and possibly 
higher-order effects in e, such as the departure of Sf of the thermal distribution functions 
from their equilibrium values. In short, it is apparent that EWB is not yet a solved problem, 
but rather one that calls for additional study. 

Undertaking this effort will be important for electroweak phenomenology. As illustrated 
here as well as in other studies {e.g., ^3^), determining the viability of EWB within a 
given electroweak model involves a detailed interplay of collider phenomenology, precision 
electroweak data, EDM searches, and a careful treatment of the dynamics of the electroweak 
phase transition. In particular, in light of the open questions pertaining to the latter, it is 
too soon to draw definitive conclusions about the implications of the next generation of 
EDM searches for the baryon asymmetry. One hopes, however, that by the time these 
searches obtain their first results, the context for their theoretical interpretation will have 
been further clarified. 
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^ An examination of Eqs. (|B1IB3|I and l|B5IB7|l of Appendix^l etc. indicates the presence of an additional, 
dynamical enhancement factor ~ ui/ JAu;'^ + F?. in the relevant integrals. 
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APPENDIX A: PROPAGATORS AT FINITE TEMPERATURE AND DENSITY 



In this section, we derive some useful properties of propagators at finite temperature and 
density, using derivations based on those for the case of finite temperature and zero density 
in Refs. Hlli. 

1. General Structure of Fermion Propagators 

We begin with the spectral function for fermions at temperature T = in the presence 
of a chemical potential /i: 

p^p{x) = Ixr [e-^(^-^^){^«(^),^/3(0)}] , (Al) 
where Z = Tr[e^^'-'^^'^'^'']. It is convenient to define the retarded and advanced propagators: 

S^(x) = ^(x°)p(x) (A2a) 
S^{x) = -e(a;°)p(x), (A2b) 

supressing spinor indices. The Fourier transforms of S^'^{x) and p{x) are related by: 

^''(^■°^)-'r^^ (A3a) 



■oo 
oo 



S\k',k) = ^| pPt2R_. (A3b) 
^ ^ 27r fco - - ze ^ ^ 

It is possible to express the momentum-space spectral function in terms of a single product 
of ipa{x) and i)p{x) instead of the anticommutator in Eq. ()A1|) . whose Fourier transform is: 

r 1 (A4) 

n 

Now insert a complete set of states between the fermion fields: 

(A5) 

+ (n| e-^(^^-^^)^^(0) \j)m^{x) \n) . 

We can rewrite the second term by switching summation labels and translating ip^ from x 
to 0: 

^^^|^-/3(H-,N)^^(0) |j) {j\i,^{x) \n) 



^ (A6) 

^gi(£„-i?,)tg-i(k„-k,).Xg-/3£,g/3MA^„+i) ^^(0) |j) (j| ^^(0) \n) , 
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which after integrating in Eq. ()A5|1 . becomes 
i5^(27r)^5(a; + E„-E,)5'(k + k„-k,)e-^(^"+-)e^'^(^"+i) (n| .^^(O) \j) (jI^^O) \n) , (A7) 

where we used the first delta function to replace Ej with En + u in the exponential e~^^^ . 
This can now be written: 

which is e~^^'^~^^ times the first term of Eq. ()A5|) . so we conclude: 

p(cu, k) = [1 + e-^("-^)] y rf^x e*M-k.x)^ ^g-/3(H-MN)^^^3,)^^(0)] . (A9) 

Similarly, we could have manipulated the first term of Eq. ()A5|) in the same way, and derived 
the companion relation: 

p{uj,k) = [1 + e'^('^-^)] J d^xe'("")^Tr[e-^(^^-^^)^/3(0)^,(x)]. (AlO) 

The Green's functions S^{k^,k) and — S'^(A;°,k) appear on the right-hand sides of 
Eqs. flA9IA10|) . giving the relations: 

5>(fc°, k) = [1 - npik'^ - /i)]p(A:°, k) (Alia) 
S<{k\ k) = -npik'^ - li)p{k\ k), (Allb) 

where npi^x) = 1/(1 + 6^). 

The various Green's functions satisfy the identities: 

S\x, y) = y) + S< {x, y) = S\x, y) + S> (x, y) (A12a) 

S\x,y) = S>{x,y) - 5«(a;,y) = S<{x,y) - S\x,y) , (A12b) 

which follow directly from the definitions in Eqs. ()8IA2|) . Thus, using Eq. ()A11|) . the time- 
and anti-time-ordered propagators can be expressed in terms of the retarded and advanced 
propagators: 

S\k\ k) = [1 - - p)\S\k\ k) + - p)S^{k\ k) (Al3a) 

S\k\ k) = -nF{k^ - ^i)S''{k\ k) - [1 - npik^ - p)\S^{k\ k) . (Al3b) 

Also note that p = - = S> - S< . 



2. Bosonic Propagators 

Similar results may be derived from scalar bosonic propagators, for which the analog to 
Eq. (0^ is: 

k) = [1 + UB^k"^ - p)]p{k°, k) (Al4a) 
G<{k^,k) =nB{k'^ - p)p{k\k) , (Al4b) 
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where the momentum-space spectral function p(/c'', k) for bosons is the Fourier transform 
of: 

p{x) = iTr{e-'^(^^-'^^)[0(x),0*(O)]} . (A15) 
The bosonic propagators also satisfy the identity p = — = G^ — G^. 



3. Tree-Level Propagators 

At tree level, the propagators S^'^ for fermions are given by: 



and G^'^ for bosons are given by: 



where = |k| + m?. Note that these propagators are independent of the temperature 
and chemical potential, which only enter in the thermal distribution functions appearing in 
the relations of the retarded and advanced propagators to the other Green's functions, for 
example, in Eq. ()A13|) . 

4. One-Loop Corrections to Massless Fermion Propagators 

Resumming the one-loop self-energy into the fermion propagator at finite temperature 
changes the pole structure of the propagator dramatically, introducing a new collective 
"hole" excitation of the plasma 0, In fact, this structure can be shown to hold even 
beyond perturbation theory Extending the results of Ref. jl^ to include dependence 
on a chemical potential, the propagator takes the form given in Eqs. (1281 - 0111) . Recall that in 
those equations Sp^h = ^p,h — i^p,h are the complex poles of the spectral function, and Zp^h 
are the corresponding residues. At leading order in the "hard thermal loop" approximation 
(see Ref. ^J), calculating the poles only to order £ ~ gT, one finds T = 0, and Zp^h{k,fi) 
and u!p^h{k,fj,), where k = |k|, depend only quadratically on fi/T, which we thus neglect in 
our analysis in the present work, where we keep only effects linear in p/T. In this limit, and 
including only a single gluon loop in the quark self-energy diagram, the poles of the spectral 
function are given by the solutions to the equation: 



^ ,0 , a.C^vrTHA k^\ . k^ + k 
= k — k — 



4:k 



T 



log 



fco - A; 



(A18) 



where Cp = 4/3 is the Casimir of the fundamental representation of SU{3). The solutions 
to this equation give the poles k^ = Ep{k), —Eh{k). The residues satisfy: 

^2 _ ^2 

Zp,h{k)= , (A19) 
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where 



rrij 



(A20) 



Calculation of the imaginary parts Vp^h of the poles, since they begin at order g^T, requires 



a resummation of hard thermal loops in self-energy diagrams jl9l l20l . |21|. We are also 
interested in their dependence on the chemical potential /i. We leave the calculation of 
these effects to a future study. 



APPENDIX B: EXPANDED SOURCE TERMS FOR QUANTUM TRANSPORT 
1. Bosons 

The CP-conserving source term for right-handed stops in Eq. ()42j] can be expanded by 
explicitly taking the imaginary part of the integrand: 



1 ^cUt _^*^^(a;)|2 



T 27r2 



OO 7,2 



k'dk 

^R^L 



(Bl) 



— (sin (f)lmh~^ + cos Re /i^) — - (cos 6Reh^ — sin 6 Im /ij) 



— (sin (plmh^ — cos Re h^) + - (cos 9Reh~j^ — sin 9 Im /i^) 



where 



A 

6 

tan^^ 
tan 6 



kr + M? 



(B2) 



v/(ri + rij)2 + (a;,.-a;R)2 

UJL + Co-R 

exp[(u;L,R ± iTL,R)/T] 
{exp[K,^,±^ri,R)/T]-l}^ 



and where are the thermal widths for the ti^R. The rates defined in Eq. can 
then be expressed: 



""^ PdA; r 1 

X 



T47r2 



(B3) 



t^RCUL l^A 

1 

+5 



[sin01m(/i^ ± /i^) — cos0Re(/;.^ =F /i^ 
[cos 9 Re{h\ =F Z?-^) - sin 9 Im(/i^ =f /ijj 
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Meanwhile, the CP-violating source given in Eq. ()48|) can be expanded: 

Sf^{x) = Ncyl\r^{liAM^fK^) ^^^'"^^ ^ 
1 



27r2 UlUJr 



(B4) 



52 
1 

A2 



[Re (1 + ra+ + n+) sin 29 + Im {n^ + n^) cos 29] 



+ — [Re (n^ — nj^) sin 20 — Im (n^ + nj^) cos 20] j- 



where ri^ ji = nBiuJi^^ ^ =t Tl^r). Our result agrees with that of Ref. js| except for a different 
relative sign in front of the cos 20 term and the overall factor of Nc- 



2. Massive Fermions 



The CP-conserving rates for Higgsino-gaugino interactions given in Eq. (jHUIl can be ex- 
panded: 



J 



where 



27r2 yuOTjUJ 



(B5) 



X 



^^^^'^ ^ii^w -k'^ + M2 cos6'^ sm2p{x 



X 



cos Re(h~ =F hi) — sin lm(h~ ± /it ) 



+ [r j^w^ - ^w^h] sin Re(/ii T /i± ) + cos Im(/ii ± /it ) 
-F^l [uJhUJ^^ - ^h'^w + - M2 cos6'/,sin2/3(x) 



X 



cos 9 Re(/ii T h±) - sin ^ Im(/ii ^ /it ) 



- [r^o;^ + T^Uh] [cos lm(/it, ^ /i+ ) + sin ^ Re(/i± ^ /it ) 



A 

5 

tan ^ 
tan0 



/i± ~ 



(B6) 



{exp[(w^^ s ± ^r^^)/T] + ly 
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The CP- violating Higgsino source in Eq. ()64p can be expressed: 



sf^{x) = 2glM^\ui{li)v{xfp 



OO 7,2 



k'dk 



X 



1 

A2 

1 

52 



w 



H 



sin 2(bRe[N± - N±] + cos 26 Im A^± + A^: 



W H 



(B7) 



sin 2^ Re[l- N± - N±] - cos 2^ Im A^± + Nt 



w 



H 



W 



H 



where -^^^ = '^b{'^hw ^■'^5 vk)- result agrees with that of Ref. ^ except for the 
sign of the cos 20 term. 



3. Chiral Fermions 



For chiral fermions, the CP-conserving chirality-changing rates in Eq. can be ex- 
panded: 



_ 1 NcVtvl 



tR 



where 



T 



X 



k^dk 



(B8) 



sin^p{Re(Aj/i+^ T A^h) - MKl ^ Kr)} + cos^pRe(/iJ^ T K^) 



T 



+ - cos2^p(Aj T Af ) Re(l - iV+ - iV+,) 

Op 



7L 7i? 



A 



hp 



sin0;,p{Re(Aj/i+^ ± A^/i+r) - Im(/i+^ ± /i+j)} - cos 0;,^ Re(/i+^ T h 
T 



hRJ 



+ 



A 



hp 



■ COS 20,p(A,^ ± A^) Re(iV+ - iV+j) 



+ {p ^ h) 



id. 



R 



L^2 



^?L = /^FK''±0,etc. 



tan Qr 



tan ( 



and where the 



n^(cj^±ir^), etc. 
^p + < 

p ' p 
- ^p 



xL,R _ p,h 



(B9) 



(BIO) 
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parameterize the hnear shifts in the thermal widths due to non-vanishing chemical potential. 
As noted at the end of Appendix ^ in a fully resummed calculation of the fermion self- 
energy, such shifts which are linear in iii/T may arise, and thus have to be included in 
our calculations, which we defer to future work. Also note that we have approximated the 
residues Z^'j^ to be purely real, which is true at the order we are working. 
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